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ABSTRACT 


The development of the free electron laser (FEL) as a source of coherent 
radiation continues toward the fulfillment of its potential for high power, high efficiency 
and short wavelength. New expenments toward fulfillment of the FEL’s potential 
present new phenomenon to be studied by theoreticians and experimentalists. Two of 
these phenomenon in short wavelength FELs are the shift in phase velocity resonance 


and the reduction in gain. 


The first part of this thesis presents an argument for the use of the FEL in a 
mantime military capacity. The advantages of the FEL over conventional kinetic 


systems and other high-energy laser systems are emphasized. 


The remainder of this thesis examines several effects of short wavelength FELs. 
In particular, Chapter IV examines the characteristics of the proposed CEBAF UV FEL 
and presents alternative designs to assist in design selection. Chapter V analyzes the 
resonant phase velocity shift and loss in gain that occurs in short wavelength FELs 
and presents a method to predict the resonant phase velocity and gain. 

Chapter VI presents a new tool for examining the optical mode within the FEL. 
By removing the input optical field from the total optical field, only the excited optical 
field amplitude and phase remain. Chapter VI presents a modification to a 


self-consistent three-dimensional simulation that will display the excited optical field. 
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l. INTRODUCTION 


The free electron laser (FEL) uses the energy of a relativistic electron beam 
to produce coherent radiation and promises high power, efficiency and 
reliability. The FEL was first proposed by John Madey in 1970 [1] and has 
since been the object of significant research and development efforts at 
universities, national laboratories and private industry around the worid. 
Theoretical understanding and technological demonstrations of FELs are 


approaching the full potential of the FEL as a source of coherent radiation. 


There are many reasons that such Significant interest is being shown in FEL 
development. Among those reasons are the potential for multi-megawatt 
average power output, very high wallplug efficiency, and optical wavelength 
tuneability over an order of magnitude. A system of such potential shows 
promise for niilitary applications as well as significant medical and scientific 
research. The military applications which might develop with the maturation of 
the FEL have been the subject of considerable study. Ground-based as well as 
space-based FELs have been a topic of interest of the Strategic Defense 
Initiative Office (SDIO) for quite some time and have been studied extensively. 
A shipboard FEL (SFEL) is discussed in Chapter Il with emphasis on the 
advantages such a system would have over conventional kinetic energy 


systems and other high-energy laser (HEL) systems. 


Chapter Ill gives a brief overview of FEL theory as a lead-in for the 
development of subsequent chapters. In Chapter IV, the ultraviolet (UV) FEL 
experiment proposed for the Continuous Electron Beam Accelerator Facility 


(CEBAF) is examined. In particular, the design of a short wavelength FEL 


which will take advantage of the CEBAF superconducting radio-frequency (SRF) 
accelerator is explored in the hope of assisting in the selection of the optimum 
design. The research presented in this thesis has contributed to the design of 
the CEBAF FEL. Features discovered here have led to the modification of the 
proposed CEBAF UV FEL. 


A particular characteristic of FELs that attempt operation at shorter 
wavelengths is that the optical mode and the electron beam tend to be roughly 
the same size. Chapter V presents new research on this charactenstic and its 
effect on the resonant phase velocity and gain of the FEL system. Simple 
equations show that as electron beam size increases relative to the optical 
mode size, the resonant phase velocity of the FEL is reduced. Additionally, a 
"universal curve” for predicting the gain of an FEL with a small optical mode is 


presented and described. 


Finally, Chapter VI examines the topic of optical mode distortion and 
presents a new method for examining its characteristics. Here, a numerical 
simulation is modified to present the stimulated optical field of an operating 
FEL. This simulation is a tool which may assist in the understanding of the 


optical field development in an FEL. 


ll. A SHIPBOARD FREE ELECTRON LASER 


A. MOTIVATION FOR A SHIPBOARD HIGH-ENERGY LASER 


There are several factors which make the development of a shipboard 
high-energy laser (HEL) system a logical and possibly essential progression in 
naval weapons development. Two major factors are the modern threats which 
face the surface ship battle groups (BGs) and the addition of anti-ballistic 
missile (ABM) and anti-satellite (ASAT) missions to the missions of sea control 


] 
and power projection. 


In the past decade, quantum advances in the performance of anti-ship 
missiles (ASMs) have led opponents of surface combatant warfare to conclude 
that surface ships (and aircraft carriers in particular) are becoming obsolete. 
These opponents surmise that the advent and proliferation of sophisticated 
long-range, high-speed, low-flying missiles make protection of BGs so difficult 
that they are virtually indefensible. Missiles fired from short range or missiles 
incorporating stealth technology will require a rapid response from a ship’s 


defensive systems. 


In an interview for the Navy Times, Commander Gregory Dundas, Deputy 
Head of the Antiair Warfare Division in the office of the Assistant Chief of Naval 
Operations for Naval Warfare, said, "Ship self-defense is our greatest challenge 
and most immediate weakness." [3] Quick Reaction Combat Capability 
(QRCC) is a high-priority program to improve surface ship defense against 
sea-skimming, antiship cruise missiles. Systems like QRCC will become 


increasingly important as force reductions take place. The same Navy Times 


article reported that Navy officials said, "In the future, deployments will be 
responsible for defending larger areas of the ocean," and Dundas added, "Sea 


control is a must. We need to do more with a smaller force of ships.” [2] 


Statements like these demonstrate the need for a more potent naval 
weapons system. With the additional threat of modern theater or tactical | 
ballistic missiles (TBMs), the argument is even stronger. An HEL, and in 
particular a free electron laser (FEL), if applied to ASM and TBM defense, 
would be a quantum improvement over the current QRCC plan, and would be 
the solution that would enable the fleet to successfully defend itself and 


continue to carry out its vital missions. 


One of the difficulties associated with defending the BG or an individual 
ship is the time lag between the detection of a threat and its destruction. The 
most significant part of this time lag is the time between the employment of a 
defensive surface-to-air missile (SAM) system and the destruction of the threat. 
In the engagement of a supersonic ASM at long range with modern SAMs this 
time lag can exceed 2 minutes. Even at short range, the time between the 
defensive system employment and threat destruction can be more than 
10 seconds. In the time it takes a SAM to transit to the ASM, the ASM may 
have moved more than 3 km closer to the ship and a second shot, if necessary, 
may not be possible. A weapon that fires lethal energy at the speed of light 
has a tremendous advantage over kinetic weapons in _ shoot-look-shoot 


sequences. 

A second difficulty associated with defending a ship or BG is that of limited 
SAM inventory, known as magazine depth. Due to probability of kill (PK) for 
each SAM, hostile electronic counter-measures, and deception, it is likely that 


multiple SAM engagements would be required to destroy each inbound ASM. If 
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faced with a large, coordinated attack with perhaps dozens of ASMs, a BG 
could quickly expend all of its available ordnance. The last ditch defense of 
even the Close In Weapons System (CIWS) is limited to only a few 
engagements, and even an ASM kill at short range may result in significant 
damage to the defenders due to debris from the ASM. Therefore, a defensive 
system which has a limitless supply of “ammunition” has significant advantage 


over present and projected future defensive missile systems. 


The Strategic Defense Initiative (SDI) began in a speech March 23, 1983 
when President Ronald Reagan called on the United States scientific 
community to develop a syste n which, ”. . . could intercept and destroy 
Strategic ballistic missiles before they reach our soil. . . ." [3] High-energy 
lasers are expected to play an important role in ballistic-missile defense (BMD) 
and FELs, as well as excimer lasers, are the most attractive candidates for this 
mission [3]. As a part of SDI, FELS were proposed as HELs that could be used 
as ASAT weapons cr as a defense against intercontinental ballistic missiles 
(ICBMs). Many organizations have worked to develop a space-based FEL 
(SBFEL) concept including Los Alamos National Laboratory, Lawrence 
Livermore National Laboratory, TRW, and the Lockheed Missile and Space 
Company. One such concept was developed by Lockheed under contract to 
the Air Force Space Division [4]. A TRW point paper suggested that the 
Capabilities and technology required of the SBFEL could be adapted to the Ship 
Missile Defense FEL (SMDFEL) concept [5]. A combination of the SBFEL and 
SMDFEL concepts could fulfill significant needs of the Navy including fleet 
defense against strategic and tactical aircraft, cruise and tactical anti-ship 


missiles and tactical ballistic missiles. 


A concept presented to the Physics department of the Naval Postgraduate 
School by Lieutenant Colonel Ed Pogue, USA, former Deputy Director of 
Directed Energy for the Strategic Defense Initiative Office, coupled a shipboard 
FEL with a high-altitude mirror relay. This concept, named Thunderball, was 
studied by Jim Bell of the Advanced Technology Group of W. J. Schafer 
Associates. [6] The high-altitude relay would provide for the additional missions 
of cueing and designation for ground-based kinetic energy weapons, contingent 
theater defense against theater ballistic missiles, tactical satellite interdiction, 


and limited or full global protection against limited strikes (GPALS) [6]. 


B. ADVANTAGES OF THE SHIPBOARD FEL 


There are many charactenstics that make the SFEL a logical choice for 
development as a naval weapon. Among these charactenistics are the 
speed-of-light delivery of lethal energy, the relatively infinite magazine, the short 
and tuneable wavelength, the high power scalability, and the very high 
efficiency. 

The speed of light is the fastest possible speed for the delivery of any 
type of lethal energy. The speed of light delivery of lethal energy by an HEL 
would greatly reduce the time lag associated with the time of flight of SAMs. In 
either an offensive or defensive role, a light-speed weapon could rapidly engage 
many targets in a very short time and, with continuous engagement, would have 
a high PK. As an example, suppose an ASM traveling at Mach 2 is detected 
and identified by a ship at a range of roughly 10 km. From the time of launch, 
a SAM traveling at Mach 2.5 would require about 7 seconds to intercept the 
incoming ASM. By this time, the ASM would be within 5 km. After the 


estimated time of intercept, an assessment of kill or miss would have to be 
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made before a re-engagement could be attempted. If the shot was evaluated 
as a miss, a second SAM would have to be fired or, if the target was at very 
short range, the CIWS would have to be engaged. Realistically, there is no 
time remaining for re-engagement and a kill by CIWS at such a close range 
would probably result in damage to the ship. In contrast, the use of a light- 
Speed weapon would require about zero transit time (roughly 0.3 ps). 
Depending on the HEL design, the laser dwell time for destruction could be as 
short as one second, and the HEL could continue to engage the target until the 
target is evaluated as destroyed. There is virtually no delay between ASM 


engagement and destruction. 


The all-electnc nature of the SFEL complements the Navy's integrated 
electnc drive (IED) concept for the all-electric ship of the future. For instance, 
an FEL with 1 MW average power operating at 25% walliplug efficiency would 
require 4 MW of prime electrical power. An IED design created by General 
Electric [7] utilizes two LM2500 gas turbines that deliver power to two 
propulsion generators. Each generator is rated at 22,187 KVA or 22.2 MW. A 
provision is also made for two 2500 KW propulsion derived ship service (PDSS) 
generating systems. Thus, this dual or twin turbine system produces a total 
electncal power of 49.4 MW and operates at 93.6% efficiency. With larger 
systems (for instance four rather than two gas turbines and associated 
equipments) the total output power of 98.8 MW is very feasible. Although this 
would require a larger ship to accommodate the larger engineenng plant, the 
technology is definitely scalable. For instance, the Spruance class destroyers 
and the more modern Ticonderoga class cruisers currently used in the Navy 
operate with marine gas turbine propulsion plants with four LM2500 gas 


turbines. These ships also use several smaller marine gas turbine electncal 


generators. The addition of superconducting technology to the IED system 
provides potential for even greater power outputs and efficiencies. Because the 
FEL is an all-electric system, it has the advantage of having a virtually limitless 
supply of “ammunition”. The FEL depends only on the amount of fuel a ship 
can carry to maintain a steady supply of electrical power not a magazine of 


“bullets” or weapon specific chemicals like a CO, laser. 


An FEL'’s optical wavelength is governed by the resonance condition [8] 


r= 2, 1+K) : a 2-1 


where A) is the undulator wavelength, K is the undulator parameter and y is the 
electron beam Lorentz factor. A broad range of wavelengths and easy 
tuneability are two factors that set FELs apart from conventional lasers that are 
tied to the natural resonance frequency of the atom or molecule. An FEL can 
be designed to operate in any wavelength band from nm to mm by selecting an 
appropriate undulator wavelength and electron beam energy. An FEL can be 
designed to take advantage of a wide range of wavelengths. Ht is possible to 
design an undulator with a range of undulator wavelengths available and it is 
fairly simple to alter the energy of the driving electron beam and thus the 
Lorentz factor. It is also possible and quite simple to adjust the undulator 
magnetic field, and thus the undulator parameter, to tune an operating FEL to a 
very specific wavelength. FELs have demonstrated operation from 240 nm 

to 18 mm and have demonstrated continuous tuneability of a single FEL over 
the operating wavelength by a factor of 10 [8]. The FEL can therefore be 
designed to adapt to new understanding of the most favorable atmospheric 


transmission band and counter-measures, and can be tuned to account for local 


atmospheric variations to obtain optimum transmission. Still, atmosphenc 
propagation at low altitudes may be the most difficult obstacle for the 


application of FELs to naval missions. 


The absorption of optical radiation by the atmosphere is dependent upon 
wavelength in a very complex manner. This absorption leads to thermal 
blooming. Fortunately, there exist certain wavelength regions with low 
atmosphenc absorption. A calculated, high-resolution, atmospheric 
transmission spectrum at 3.8 ym over 10 km sea level path at 0°C and 76% 
relative humidity shows transmittance of better than 95% [9]. Transmission 
at 3.8 um would minimize the ef’ects of thermal blooming. Though aerosol! and 
turbulence effects become dominant and lead to large variations in performance 
at the shorter wavelengths, the 3.8 num wavelength offers the best performance. 
At longer ranges performance can be improved by using multiple-pulse 
propagation as well as active or adaptive optics techniques for turbulence 


compensation and for thermal blooming correction. [10] 


The potential for high power in an FEL has been demonstrated by several 
experiments, but the most important characteristic of a weapons grade FEL, 
high-average power, has yet to be demonstrated over long periods of time. To 
date, average powers of only a few watts have been observed in Several 
expenments [8]. 

The typical FEL optical output consists of a macropulse made up of a train 
of short, high-power micropulses. This macropulse and train of micropulses are 
a direct result of the pulsed electron bunches input to the FEL undulator from 
the accelerator. The average power of the macropulse is determined by the 
peak pulse power, the pulse length and the micropulse repetition frequency. 


The average power of the FEL is determined by the average macropulse power 
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and the macropulse repetition frequency. While high peak power and high 
micropulse repetition frequency have been demonstrated, the technology for a 
continuous train of micropulses is still being investigated. The Continuous 
Electron Beam Accelerator Facility (CEBAF) UV and IR FEL experiments will 


demonstrate continuous electron micropulses and high average power [11]. 


High peak power has been demonstrated in several FEL expenments. 
The ELF Il microwave FEL at Lawrence Livermore National Laboratory (LLNL) 
was able to produce 1 GW peak powers over 50 ns pulses [12]. The Los 
Alamos FEL produced peak intercavity power of 2 GW with peak output power 
of 40 MW over 16 picosecond micropulses at 10 pm wavelength [13]. Neither 
of these facilities, however, have been able to demonstrate the continuous 


micropulse repetition frequency necessary for high average powers. 


To date the highest average power achieved is 10 W produced by the 
superconducting accelerator (SCA) FEL experiment at Stanford [8]. The SCA 
expenment used 5 ps duration micropulses at 12 MHz with 100 ps duration 
macropulses at 20 Hz. This average power was achieved with peak power of 
only 1MW. The proposed CEBAF UV FEL may demonstrate average power at 
least 2 orders of magnitude higher [11]. With a continuous train of 1 ps pulses 
at 7.5 MHz and peak power of 480 MW, the CEBAF UV FEL should 
demonstrate average powers near 4 kW. Continued research into 


superconducting RF accelerators will probably produce even greater results. 


The wallplug efficiency of an FEL will be defined as the average optical 
power output of the FEL divided by the prime electrical power input. The 
extraction efficiency of an FEL is the fraction of the electron beam power that is 
converted into optical power in one pass through the undulator. While the 


output power of the laser depends upon the extraction efficiency of the FEL 
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undulator, the wallplug efficiency is more dependent upon superconducting 
technology and electron beam recirculation. Both TRW and CEBAF estimate 
that wallplug efficiency could be as high as 40%. The development of 
superconducting radio frequency (SRF) accelerators has proven to be a huge 
technological achievement. SRF accelerators provide the capability to 
accelerate higher peak and average electron currents which in turn provide high 
gain and high-average optical power. Additionally, the high Q of the cavities 


means there are few losses as the electrons are accelerated by the RF field. 


The development of energy recovery by RF accelerators is very exciting. 
By recirculating the electron bee n through a decelerating mechanism after the 
beam has passed through the undulator, much of the kinetic energy of the 
electrons is converted back into RF energy for the accelerating cavities. Much 
of the energy which remains in the electron beam after passing through the 
undulator is not lost to a beam dump. The Los Alamos National Laboratory 
(LANL) FEL used a resonant bridge coupler to pass RF energy from a set of 
decelerating cavities back to the accelerating cavities [14]. The accelerated 
electron beam was passed through a FEL undulator where 0.7% of the energy 
was converted into optical power. The electron beam was then passed through 
a series of decelerating cavities where 68% of its energy was converted back 
into RF power. A similar system in which the electron beam is introduced into 
the original accelerating cavities with an approximate 180 phase shift may be 


even more efficient. 


The ability to recover energy from the electron beam is dependent upon 
the electron beam characteristics. These characteristics are in turn dependent 
upon the accelerator injector, the FEL undulator and the electron beam 


transport system. If more of the energy is extracted from the electron beam in 


p 


the undulator, it becomes more difficult to recover RF power from it. At LANL, 
lasing with about 1% extraction efficiency did not degrade the energy recovery. 
Overall, extraction efficiency of more than 5% with energy recovery of 70% in 
SRF accelerators should lead to high average power laser operations with 


greater than 40% wallplug efficiencies. 


C. THEORETICAL FEL RANGE LIMITS 


In order to establish the range limitations of a shipboard FEL, assumptions 
were made based on the best available information. In the case of the SFEL, 
the primary assumptions involved the required flux on target and especially on 
atmospheric propagation of high-power lasers. The information available which 
describes the required flux on target varies from 1 kJ/cm* to 100 kJ/cm*. 
Assuming an energy flux of 10 kJ/cm*, whch corresponds to a power density on 


target of 10 kW/cm* over 1 second, accounts for some target hardening. 


In combining the SMDFEL and Thunderball concepts into a single SFEL 
concept, it becomes apparent that there are two distinct modes or mission 
profiles. The first is direct defense against an incoming threat and the second 
is long-range intercept of a threat through a high altitude relay. In the direct 
defense mode, the SFEL engages the inbound aircraft or missile with no relay 
and the light must transit to the target near the surface. For threats at very low 


altitudes, the engagement range is limited to the horizon distance [15] 
Range (km) = 3.71xVHeight(m) __, 


where the height is the sum of the director and the target heights. For instance, 
a SFEL system with beam director height of 20 m engaging an incoming target 


at 6 m is limited to a horizon range of about 19 km. This distance will be 
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extended for targets at higher altitudes, but is limited by beam director height 
for very low altitude missiles. In this direct defense mode, the optical beam 
must propagate near the sea surface where the propagation effects of 


attenuation and thermal blooming are most severe. 


A spot size (diameter) on target of 10 cm balances the requirement to 
keep power density high while staying within limits of pointing and tracking 
accuracy at long ranges to achieve the necessary dwell time. The beam director 
diameter required to focus to a 10 cm spot at 20 km is a direct result of 


diffraction limited beam so that 
D=RANn/A = Rr = (20 km)(3.8x10% m)/.05 m = 1.5 m. 


For higher altitude target, range is limited by ability to focus the diffraction 
limited beam to a 10 cm spot with a given beam director size. Given a 


reasonable director size of 3 m the maximum range is determined by [3] 
R* = AD*/nr° 
R=rD/X=40km_. 2-2 


The average laser output power required to deliver 10 kJ/cm* over a 10 cm 


spot at 40 km comes from the equation for fluence [3} 
F=PIDnR?r*  , 
So that 
P=nFR°7/tD* = 807kW_ i. 2-3 


Assuming a maximum of 25% transmission, including beam blowup losses, the 


required FEL output is about 3.2 MW. 
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There are many options available for a high-altitude relay system to 
extend ranges against high-altitude targets. Placing a relay mirror at 80,000 ft 
(25 km) avoids the most significant portions of the atmosphere and significantly 
improves laser propagation and the lethal range [6]. A larger output mirror at 
this high altitude would allow focusing over much greater distances. In this 
case, maximum lethal range is limited by laser output power and output mirror 
diameter. If the maximum power output is limited to 10 MW then, neglecting 
mirror and transmission losses, a maximum range for high altitude intercept can 
be calculated. Assuming a 6 m relay mirror and a 3 m beam director, with the 


same fluence required on target, the maximum range becomes [3] 
R?=PID*inFr? 
R=NVPtinFD/X=300km_. 2-4 
To get the laser power up to the high altitude relay requires the ability to focus 
the ship’s output to a size smaller than the input mirror of the relay. Because 
the laser output beam is nearly diffraction-limited, the minimum spot size is 


determined by range, wavelength and director diameter. A 3.8 ym laser with 


a 3m director can easily focus on a 3 m input mirror for the relay at 25 km. 


D. DISCUSSION 

The need for a very fast and effective defensive system for ships and 
battle groups has been amply demonstrated. Additionally, the Navy’s 
Capabilities must grow to meet the needs of the new ASAT and ABM missions. 


It is therefore essential that long term research be devoted to develop FELs and 
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demonstrate their full potential. While there are still many technological 
roadblocks to the deployment of the FEL as a weapons system, the proven 
advantages, the scalability, and the potential efficiency make the FEL a logical 


choice for development. 
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lil, FREE ELECTRON LASER THEORY 


A. BASIC FREE ELECTRON LASER PHYSICS 


An FEL consists of two major components. An electron accelerator 
provides a stream of relativistic electron bunches to an undulator which has a 
periodic magnetic field to “wiggle” the electrons as they pass. The penodic 
acceleration of the electrons causes them to radiate in a forward cone. Some 
of this spontaneous radiation may be saved in a laser resonator formed by 
placing two curved mirrors beyond the ends of the undulator. As the radiation 
builds, the coupling of the optical field in the undulator with the wiggling electron 
beam leads to stimulated emission and coherent radiation. Though this is a 
quantum description of the electron/optical field interaction, FELs can be 
described as classical devices and may be described with electro-magnetic 


theory [16]. 


The radio-frequency (RF) electron accelerator yields picosecond long 
bunches of electrons with peak current on the order of one to hundreds of 
amps. The electron bunches achieve typical energies, (y- 1)mc*, of tens of 
MeVs up to hundreds of MeVs giving Lorentz factors on the order of y= 10%. 
The radius of the electron beam is generally on the order of millimeters and 
yields an electron density on the order of p=10'*-10'%cm~. The beam 
quality, described in terms of the energy spread and emittance, can have a 
large effect on FEL gain. Energy spread is the spread in the electron velocities 
described by Ay/y. Emittance e = 7 6 where 7 is the rms initial position spread of 


the electrons and 6 is the rms initial angular spread. While the position and 
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angular spreads of the electrons can be changed in the beam transport system 
of the accelerator, the emittance remains constant [17]. The normalized 
emittance, c, =ye, is useful in Comparing the beam quality of accelerator at 
different energies. Typically, each electron bunch has a Spread in electron 
energies of a few tenths of a percent and a normalized emittance on the order 


of a few to a hundred mm-mrad. 


The undulator consists of a penodic magnetic field with linear or helical 
polarization. The magnetic field may be generated by permanent magnets, 
current-carrying coils, or a combination of each. Each undulator period, Ag, is 
typically a few centimeters lo1g with perhaps N=100 penods giving the 
undulator a length on the order of a few meters. The strength of the undulator 
may be described by the undulator parameter K = eBi)/2nmc* where e=le! is 
the charge magnitude of an electron, m is the mass of an electron, and c is the 
speed of light. The rms magnetic field strength B = B/V2 is typically a few 


kilogauss so that K = 1. 


The electron beam transverse dimension and angular spread contribute to 
the line width of the FEL. The spread in the line width is minimized when the 
contnbution of the radial spread is equal to the contribution of the angular 
spread. The equal contributions leads to the optimum condition for focusing the 
electron beam as [17] Kk)’, = 76, and leads to a matching beam radius of 
ro =Ve,aV2nK. This matching condition gives smooth propagation of the 
electron beam along the undulator axis with minimum spread in the line width. 
As the number of undulator periods increases and the optical wavelength 


decreases the FEL becomes more sensitive to emittance. 


We 


Figure 3-1 shows the microscopic features of the classical approach to FEL 
physics. At the top of Figure 3-1, electrons enter the undulator and encounter 
the periodic magnetic field where they begin to oscillate. As the electrons 
oscillate they emit photons in a forward cone. The middle of Figure 3-1 shows 
one period of the undulator expanded. The radiation travels along the undulator 
over one period of the undulator magnetic field as the electron follows a 
sinusoidal path over the same undulator period. The bottom of Figure 3-1 
shows that the effect of the optical field on the electron is dependent on the 
relative phase between the electron and the optical field. When the electron 
has a velocity component in the z direction (along the undulator axis) it may 
experience a retarding or an accelerating force that causes the electron to give 
up or gain energy. Because energy conservation must apply, the energy given 
up by an electron is stored in the optical wave leading to amplification of the 
Optical field. The energy can be absorbed by an electron depending on the 
phase between the electron and the optical wave. Then, the optical field gives 
energy to the electron resulting in a loss in the optical field and acceleration of 
the electron. 

In an initially random beam, spread over many optical wavelengths, both 
absorption and loss occur in roughly equal amounts. At a_ particular 
wavelength, the energy transfer to the optical field can be made to dominate 


giving rise to net gain [16]. 
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Figure 3-1: Major components of the FEL illustrating the fundamental 
physics of the laser. 
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B. ELECTRON DYNAMICS AND THE PENDULUM EQUATION 


The forces acting on the electron in the undulator are governed by the 


relativistic Lorentz force 


d & pe ceria -2 2 
SOB) . © (B+ Bx (8, +B,), += BE) y2=1B-B . 3 


where V=fec is the velocity of the electron, E, = E[cosy,-siny,0] and 
B, = E[siny,cosy,0] are the electric and magnetic fields of the optical wave and 
B,, = B[cos(koz),sin(koz),0] is the undulator magnetic field for a helically polarized 
undulator. Equations for the trajectories of the electrons in the undulator are 
derived by inserting the equations for the electric and magnetic fields into 
equation 3-1. Because there are five component equations in equation 3-1 and 
Only four unknowns [X(t), y(t)], one of the component equations can be ignored. 
When v=c, B,=1 and 2(t) = Boct +.... and the electrons travel mainly along 
the axis of the undulator. In this case, the transverse optical fields in the first 
equation nearly cancel so that the transverse optical force is much less than the 
undulator field force, IE,|(1-B,)<«<B,B [16]. When this is the case, the 
transverse force on the electrons is primanly determined by the undulator field 


and the transverse component of equation 3-1 becomes 
Tug) = -© 8, B[-sin(koz),cos(koz),0] 3-2 
ot mo‘ 7 Cad mia 


The transverse velocity of the electrons is found by integrating equation 3-2 by 


inspection. The result is 
K , 
B= Fg teos( Noa) SIN care ; 3-3 


where K = eBA,/2nmc* is again the undulator parameter. The constants of 


integration have been eliminated by assuming perfect beam _ injection. 
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Integrating once more yields the trajectories of the electrons in the undulator as 
K 
Xx, = = -sin{aot),cos( cot) 0 ; 3-4 


where @) = Koc. Because yy 1 and K = 1, the transverse oscillations are small 


(AX /Ag = K/2xy) compared to the undulator wavelength. 


Because the magnetic field increases off-axis, a second, slower oscillation 
called betatron motion is superimposed on the oscillation of the electrons in one 
period of the undulator. The wavenumber of this motion is kg = Kko/y. This 
motion is characteristic of the natural focusing of the undulator field [16]. The 
number of betatron oscillations is Ng = kgl/2n = KL/yAg = KN/y. If Ne is small, 
there is no significant focusing in the FEL. [16] 

The laser gain can be calculated from the second equation of equation 3-1 
as the change in energy dy/ct = y. Substituting the optical electric field E, into 


equation 3-1 gives 





s aie e : : 
y= —(B, ae 3-5 
because E, =0. Now, 
Y = eEK cos(C + 3-6 
: a (C+) , 


where C= C(t) = y+ Koz(t) = (k + Ko)z(t) - wt becomes the electron phase. At 
t=0, ((0) =p =(kK+ko)Z2. Because k >» ko, Co = k2o = 2nz/A. The electron 
phase relates the electron z position to the optical wavelength. 

To find the evolution of the electrons in the presence of an optical field, the 
Lorentz factor is written in terms of the electron phase using the fact that 


y* = 1-B* = 1-B?-B? and B, = K/y. When y is large, B, = 1- (1+K)/2¥*. Taking 
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derivatives of the electron phase, 
( = (k+k))z —@= koBpe—KUI-Bejo * S7 
and 


t = (k+k)z = roll) : 


14+K2 


a 


Because w(1+K*)/2;7 = Wp = ko¢ equation 3-8 can be written as 


me cos(t+p) 3-8 








C= ke 








.. EKk 
C — eNO aie) 
ym 


The dimensionless time can now be defined as t= ct/L so that t= 0-51 as 
the electrons pass down ine undulator. Now, o(..)/ct = (L/c) d(. dt = (..), and 
d*(..)/at® = (L/c)*0*(..)/dt? = (. °) so that the pendulum equation can be written as 


Liao 





C(t (t) = cos(C+p) = lalcos(C+o) 3-9 


where lal=4nNeEKLiygmc* is the dimensionless optical field strength. 
Equation 3-9 is the pendulum equation which describes the evolution of the 
electrons through the undulator. In equation 3-9, if -3x/2 < (C+) < w2 the 
electrons will lose energy to the optical field. Similarly, if m2 < (C+) < 32/2 the 
electrons will gain energy from the optical field. 

If the electron phase velocity is defined as v = a and the equation for C 
above is used, in dimensionless form v(t) = L[(k+ko)B,-k]. The electron phase 
velocity is a measure of the resonance between the electron beam, the 
undulator, and the optical field. If v(t) =0, then the FEL is at resonance and 


=(k+Ko)B, and k=koB,/1-B,. Because B, = 1-(1+K%)h* =1 for y>>1, then 
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k = 2y°k)/(1+K*) and the resonant wavelength is 


2 
= 3-10 


Equation 3-10 shows that the resonant operating wavelength of an FEL can be 
altered by changing the electron beam energy (+1)mc*%, the undulator 


wavelength A», or the undulator parameter K. 


C. THE SELF-CONSISTENT WAVE EQUATION 


Maxwell’s wave equation can be used to describe the optical evolution in 
an FEL. Spontaneous emission in an FEL oscillator grows to form a classical 
wave with a bandwidth comparable to the inverse of the number of undulator 
periods, 1/N [16]. This means that the wave has some degree of coherence 
even after one pass. The narrow bandwidth of the laser allows the assumption 
that the optical wave varies slowly in space over one optical wavelength 
(E’<<kE, o’«ko). This is the slowly-varying amplitude and phase approximation 
(SVAP). This can be viewed as a carrier wave of a single frequency which is 
modulated by a complex wave envelope that is slowly-varying in amplitude and 


phase over many optical wavelengths. 


Neglecting transverse effects of the optical field, it can be assumed that a 
circularly-polarized plane-wave which is present in the undulator has a vector 


potential of the form 
A = AD siny, cosy, 0), 3-11 


where E(t) is the optical electric field, k is the optical wavenumber and 


y = kz-wt+. Using the vector potential in Maxwell’s wave equation [16] 
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v ~ 42 )a--4 —h C 3-12 


where J, is the transverse current from the transverse motion of the electron 
beam in the Coulomb gauge [16]. The SVAP is used to eliminate all spatial 


terms with two derivatives or higher to get 
2 1 & |p 2 
| ce? 7H? 


The fast rotation of the sine and cosine functions of y can be eliminated by 








cosy, oF AT, 3-13 


defining two orthogonal unit vectors, &, ={cosy, —siny, 0] and &, ={siny, cosy, 0]. 
By projecting equation 3-13 onto ¢, and ¢2, Maxwell's wave equation becomes 


two first-order scaler equations, 





2 1 2 2 dE, 
0 - 45) C Ot =-“EY, €, a 3-14 
and 
1 0° 2E, 9 4n 
at |r a — 


The single-particle current is J, = —ecB53(% - 7) where 7 is the position of 
the i” electron and 5°(% — 7) is the three-dimensional Dirac delta-function [18]. 
A volume element, dV, which is much smaller than the coherence volume but 
much larger than the optical wavelength, is selected. By substituting the 
single-particle transverse current into equations 3-14 and 3-15 and summing 
over all the particles and averaging both equations at a fixed time over dV [16], 


equations 3-14 and 3-15 become, 


E 
aoe = piri <cos(C+9)>  , 3-16 


c at 
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and 


E+ a: omen csin(t 2 3-17 


where p is the electron particle density and <.> represents the ensemble 
average over all the electrons. Equations 3-16 and 3-17 are simplified by 
introducing the dimensionless current j = 8N(enxKL)*p/ymc* [16]. With the 
dimensionless current and the dimensionless field strength, equations 3-16 


and 3-17 simplify to 
lal = —j<cos(C + 9)>, 6 — int +o)>  , 3-18 


or in phasor form 
a =-j<e'S>  , 3-19 


where a = lale’® is the complex dimensionless field. Equations 3-18 and 3-19 
show that the bunching of electrons in phase at C=x drives the optical 
amplitude and leads to gain. While the optical phase 9(t) is driven when the 
electrons are bunched at C=72. Growth of the optical wave increases with / 
and is dependent upon the electron distribution <..>. Usually the electron 
bunching is not perfect and drives both the optical amplitude and phase. 

In reality, the evolution of the optical wave in the FEL can be significantly 
affected by transverse effects such as diffraction. To include transverse effects, 
the three-dimensional FEL wave equations are derived similar to equations 3-18 


and 3-19 from the parabolic wave equation [19] and becomes 
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ig2, 0 ei : 
-i9 + a a(X.th=—<jet> . 3-20 


Equation 3-20 can be used to perform numenical integration to simulate an FEL 


including the effects of beam size, beam quality, and diffraction. 


D. FEL DIMENSIONLESS PARAMETERS 


Dimensionless parameters are useful in discussing and comparing the 
attributes of a particular FEL. Several FEL dimensionless parameters have 
been defined including the undulator parameter, K, the electron phase and 
phase velocity, ¢ and v, and the optical field strength, !a!. There are several 


more FEL dimensionless parameters which are useful. 


When discussing the sizes of the optical mode and the electron beam, it is 
useful to compare their dimensionless parameters. All transverse dimensions in 
the FEL are normalized by vz/L’. The dimensionless electron beam radius is 
O, =1Vn/Li. 

The Rayleigh length, Zo, is a measure of the optical beam diffraction and is 
determined by the optical cavity configuration including the mirror curvature and 
separation. The optical mode waist, Wo, is related to Z) by xWé =AZp>. The 
dimensionless Rayleigh length is z)>=Z,/L [17]. When the transverse mode 
radius is normalized wy = WoVW/LA = VZ>\/n x Vr/Li, SO that Wo = VZo. 

The slippage distance, As, is the distance the light passes over the 
electron bunch as the electrons travel down the undulator. Comparing the 
difference between the distance traveled by an electron and one wavelength of 
light over one pass down the undulator shows that As = Av At = (c-B,c) L/c = 


(1-B,)L = (1-B,)N%>. Because (1-B,) = (1+K?)/2y and = A_(1+K")/2y*,_ then 
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As = NA is the slippage distance. It is the characteristic length over which the 
electron and light can exchange energy during one pass. If the electron pulse 
length is much larger than the slippage distance, each part of the optical pulse 


experiences gain proportional to the local electron density. 
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IV. THE CEBAF ULTRAVIOLET FEL DESIGN 


A. BACKGROUND 


The Continuous Electron Beam Accelerator Facility (CEBAF) under 
construction in Newport News, Virginia is designed to serve as a world center 
for nuclear physics with the first operation of the full accelerator scheduled to 
begin in 1994 [11]. The Southeast Universities Research Association (SURA) 
manages CEBAF under the Department of Energy (DOE). For nuclear physics 
experiments, the electron accelerator will use superconducting radio-frequency 
(SRF) cavities to accelerate electron bunches of approximately 0.3 pC charge 
at a repetition rate of roughly 1.5 GHz. CEBAF will have the capability to 
deliver simultaneous electron beams of energies 0.8 GeV to 4 GeV to three 
separate experimentation halls. Parallel to these electron beams, CEBAF 
proposes to provide accelerated electron beams for use in an infrared (IR) and 
an ultraviolet (UV) FEL. These FEL experiments provide an excellent step in 
the technological development of high-efficiency, high-power FELs. To assist in 
the selection of a design for the CEBAF UV FEL, reference [20] was prepared 


for the 14” International Free Electron Laser Conference. 


CEBAF provides an superb power source for a UV FEL oscillator. [11] 
The electron beam energy is (y-1)mc? = 400 MeV where m is the electron mass, 
c is the speed of light, and y= 748 is the Lorentz factor. The continuous train of 
electron micropulses is 1 ps long with peak current 120 A, energy spread 


Ay/y = 0.002, and normalized emittance ce, = ye = 15 mm-mrad. 
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The design of any short wavelength FEL requires addressing problems 
that originate primarily from limited electron beam quality. The specific features 
of the CEBAF UV FEL are typical of other short wavelength FELs. The optical 
mode waist area is proportional to the optical wavelength, A. When A is small, 
the optical mode area tends to be small, and can be smaller than the electron 
beam. The electrons outside the optical mode do not participate in the gain 
process and are wasted. The electron beam is limited in size because of finite 
emittance. This chapter examines and evaluates the performance of the 
proposed CEBAF UV FEL with recommended deviations from the ususal design 


criteria [11] in order to optimize gain [20]. 


B. THE CEBAF UV FEL DESCRIPTION 


The basic undulator for the CEBAF UV FEL is linearly-polarized with a 
wavelength of A,~=6 cm over N=50 periods and length L=Na,=3 m. 
Electromagnetic coils provide the peak undulator field of B = 4.4 kG which gives 
an undulator parameter of K = eBA)/2V2nmc* = 1.76 where the rms field is 
B=BN2. The matched electron beam size is r, = (€,A)/V2nK)' = 0.03 cm, 
which minimizes the effect of emittance on beam quality by equalizing the 
contnbutions from the radial and angular spreads. The number of betatron 
Oscillations along the undulator length is only NK/y= 0.1, because of the large 


Lorentz factor. 


The resonant optical wavelength is 4 = 4)(1+K?)/2¥° = 2000 A. Because of 
the short operating wavelength, the small optical mode _ waist, 
W, = VZoA/n = 0.03 cm, is approximately equal to the electron beam radius. 
Gain is reduced because electrons near the edges of the beam see weak 


optical fields for bunching. 
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The electron beam density is p =7x10'* cm~ giving a dimensionless 
current density of / = 8N[xeKL]*p/mc* = 2.7. There are two correction factors 
which can be included in j to account for other effects. The filling factor 
F = (r,/Wo)*(1+L°/12Z¢) accounts for the relative size between the electron 
beam and the optical mode as it diffracts along the length of the undulator. 
Because r, = Wo, then, for Rayleigh length Z)=L/2 the filling factor is large. 
The Bessel function factors, JJ = Jo(E)-J,(E), where & = K*/2(1+K?), expresses 
the reduced coupling caused by fast penodic z-motion in each undulator period. 
If € = Az/A bunching is reduced, and thus, gain is reduced [16]. The theoretical 
single-mode gain is G = 0.135/j(JJ)*F = 33%, but does not include the effects of 


self-consistent gain, beam quality, beam size, or optical diffraction. [16] 


A self-consistent gain calculation includes the changing amplitude and 
phase of the optical field as the field grows. The field causes bunching of the 
electrons and the bunching of the electrons leads to a greater increase of the 
field. When the FEL interaction is calculated self-consistently using the integral 
equation [16], the gain increases to G=36%. The length of the electron 
micropulse is /, = 0.015 cm and is 15 slippage distances in length. When the 
electron pulse length is comparable to the slippage distance each part of the 
optical pulse will experience gain proportional to the local density of the electron 
pulse. The resulting effects are known as short-pulse effects [16]. Because the 
electron pulse length is much larger than the slippage distance, there are no 
short pulse effects as is typical of short wavelength FELs. 

The electron phase velocity is defined as v = L[(k+ko)B,—-k] = 4nNAy/y, and 
measures the degree of resonance between the electron beam, the undulator, 
and the optical field. An electron beam energy spread causes a spread in 


phase velocities, og = 4nNAyy = 1.3 with Ayy= 0.002 for the CEBAF UV FEL. 
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The angular and radial spreads of the matched electron beam give equal, but 
small, contributions to the phase velocity distnbution of o,=0.3. When og, or 
Og = x, the electron beam is randomly spread over half of an optical wavelength 
during a pass through the undulator making bunching difficult and decreasing 
FEL gain. In the CEBAF UV case, the energy spread, with og = 1.3 distributed 
in a Gaussian random spread, causes the self-consistent gain calculated with 
the integral equation [16] to decrease from 36% to G=27%. The effect of 


emittance is negligible. 


C. OPTIMIZING THE CEBAF UV FEL DESIGN 


In order to maximize gain with the basic CEBAF UV FEL design, the 
effects of dimensionless Rayleigh length, z), mode waist position, t,, and initial 
phase velocity, vo, on the peak gain was examined. In the FEL interaction, the 
phase velocity and optical wavelength are not fixed and are free to evolve. 
Spontaneous emission begins on resonance at vp =0 and gain grows aS vo 
shifts to satisfy the resonance condition. A gain spectrum is obtained by 
calculating the gain over a wide range of phase velocities with a self-consistent 
numerical simulation. The peak gain can then be selected from the gain 
spectrum. Gain spectra, G(v,), were calculated over a large range of z> and t, 
to determine their optimum values and the corresponding optimum gain. The 
gain spectrum is determined by self-consistent numerical simulations which 
include the effects of diffraction, beam size and beam quality. When the effects 
of optical diffraction are included in the FEL interaction, the phase velocity for 


maximum gain is altered by the Rayleigh length [21]. In the CEBAF UV FEL, 
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the large electron beam size causes an additional change in the phase velocity. 
A detailed examination of the effects of electron beam size and diffraction is 


presented in Chapters V and VI. 


As inputs to the numerical solution for the gain spectrum, dimensionless 
parameters which correspond to the CEBAF UV FEL parameters are used. 
The peak current density is / = 2.7 at the center of the beam, the electron beam 
radius is o, = 1,¥n/L2=0.78 with a parabolic shape, and energy spread is 
described by og = 1.3 distributed as a Gaussian. The initially weak optical field 
iS aj9=1 with dimensionless optical mode radius, w?(t) = 1+(t-t,,)*/z¢, where 
t=ct/L is the dimensionless time along the undulator and t, is the position of 
the optical mode waist along the undulator. The dimensionless Rayleigh length 
is Z7=nWé/LA=Z)/L. The transverse dimensionless variables, r and w(t), 


are both normalized to VLA/x. 


For a Gaussian optical mode, a Rayleigh length of Zz) = 1/¥12 minimizes 
the optical mode volume averaged over the length of the undulator [16]. For a 
low-gain FEL with a filament electron beam, the value zy = 1/V12 = 0.3 should 
optimize the gain. An electron positioned at the optical mode waist is at 
f =Wo=VZo. Figure 4-1 shows the graph of gain versus Rayleigh length over a 
large range of values. The initial sharp increase in gain results from the fact 
that as the Rayleigh length increases the optical mode expands so that more of 
the electrons are able to participate in the gain process. The nominal value of 
Zo = 1/412, which minimized optical mode volume, is not optimum for the 
CEBAF UV FEL because of the large electron beam size. A larger value of 
Z,)=0.5 increases the optical mode waist so that more electrons participate in 


the bunching process. When the optical mode radius and electron beam radius 
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are approximately equal, the curve in Figure 4-1 begins to flatten and peaks at 
the value z)=0.5. AS Zp») continues to increase and expand the optical mode 


size, less of the optical mode sees electrons and gain is slowly reduced. 


Gain 


Figure 4-1: Gain versus Rayleigh length for the CEBAF UV FEL 
parameters. Gain peaks at Z,. = 0.5. 


Figure 4-2 is the graph of gain versus optical mode waist position, t,,, with 
the CEBAF UV FEL parameters and z, = 0.5, as determined above. The curve 
is symmetric about its peak value of t, =0.45. If the optical mode were 
positioned at the beginning of the undulator, light would quickly diffract and 
reduce the mode amplitude near the end of the undulator. The low optical 
amplitude would reduce coupling with the bunched electrons. If the optical 
mode waist were positioned near the end of the undulator, the larger optical 
mode in the first half of the undulator would provide less optical amplitude to 
bunch the electrons. In this case, the optimum mode waist position is to focus 
the light just prior to the center of the undulator to increase the electron 
bunching. Toward the end of the undulator, mode distortion helps to keep the 


optical mode focused toward the electron beam and gain is maximized. 
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Gain 
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Figure 4-2: Gain versus optical mode waist position for CEBAF UV 
FEL parameters. Gain peaks at ty = 0.45. 


Figure 4-3 shows the gain spectrum, Gj(v,), for the CEBAF UV FEL 
determined by self-consistent numencal simulation including diffraction, beam 
size, and beam quality. The optimum resonator design described by a Rayleigh 
length of z>=0.5, and waist position t, =0.45 is used in Figure 4-3. The 


maximum gain obtained is G = 20% at a phase velocity of vo = 4.3. 


Figure 4-4 shows the result of a numerical simulation of the CEBAF UV 
FEL radiation wavefront as t goes from 0 to 1 along the N=50 period 
undulator. The optical wavefront has an initially weak field amplitude a, =1 at 
the mode center with phase curvature that gives Rayleigh length 2z,=0.5 
focused at t, = 0.45 in the absence of the FEL interaction. The electron beam 
is parabolic in shape with radius o, =0.78, phase velocity vo = 4.3, current 


density j/ = 2.7 at the center of the beam, and energy spread og = 1.3. 


The evolution of the optical mode amplitude, la(x,t)l (upper-left), is shown 
as an intensity/contour plot with the x and y axes scaled to VLA/x. The 


maximum field is white and zero field is black with one contour. The greyscale 
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Figure 4-3: Gain spectrum for the CEBAF UV FEL using optimized 

parameters. Peak gain occurs at vo = 4.2. 
for each plot is shown at the nght. The mode focuses just before the middle of 
the undulator, t=, =0.45. At the end of the undulator the mode expands 
significantly in the absence of the FEL interaction, but in the plot of la(x,y)! 
(upper-center) the mode is slightly focused in toward the electron beam. The 
evolution plot of the bunching current, o(x,t) « j(r)<cos(C+)> (middle-left), 
shows the electron bunching developing along the length of the undulator. The 
final bunching o(x,y) (middle-center) is smaller than the optical mode because 
bunching is maximum where both the light intensity and beam current are 
maximum. This tends to slightly focus the wavefront even though the laser and 
electron beams are about the same size. The electron phase velocity 
distribution, f(v,t) (lower-left), shows only small distortion in weak optical fields. 


The final electron phase-space distribution in (¢,v) (lower-center) shows the 
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Figure 4-4: Numerical simulation of the CEBAF UV FEL radiation 
wavefront in weak fields. 
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slight electron bunching at the end of the undulator in the weak optical field. 
Note that some electrons remain at the initial phase velocity, vp = 4.3, because 
they are at the edge of the optical mode and see only low intensity light. The 
plot of the phase evolution at the center of the mode, 9(0,t) (on the right), 
shows the approximate linear decrease in phase « -t/z, descnbed in (2). The 
power evolution, P(t), and gain evolution, In(1+G(t)) (on the right), show the 
development of the final gain, G=20%. There is little gain in the beginning of 


the undulator while bunching develops. 


Figure 4-5 shows the result of a wavefront simulation for the same FEL as 
in Figure 4-4, but with an initially strong optical field amplitude of a, = 40. The 
phase velocity is increased to vo = 6 for more optimal gain in strong fields [16]. 
The gain evolution, In(1+G(t)), shows the gain developing early along the 
undulator due to the larger field strength. The same effect is seen in the plot of 
the bunching current o(x,t). The electron phase velocity distribution, f (v,7), 
Shows the wider spread in strong fields, while the phase-space plot, (Uv), 
shows overbunching. The bunching current evolution o(x,t) and final bunching 
distribution o(x,y) show a “hole” in the middle of the beam resulting from the 
overbunching in strong fields. Here, the electrons are actually taking energy 
back from the center of the optical mode at saturation. The small angular 
spread from the high-quality matched electron beam, 6 = e/r, = 5x10-°, does not 
allow mixing of the inner and outer parts of the beam. The transverse drift of 
an electron is only one-tenth of the beam radius over the last fifth of the 
undulator length. As the optical mode and electron distribution continue to 
evolve over many passes in stronger fields, the trapped-particle instability may 


occur and cause a hole to develop at the center of the optical mode. 
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Figure 4-5: Numerical simulation of the CEBAF UV FEL radiation 
wavefront in strong fields. 
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D. ALTERNATIVE DESIGNS FOR THE CEBAF UV FEL 


The basic CEBAF UV FEL design appears attractive because of its good 
performance at short wavelengths using the conventional FEL interaction. The 
N = 50 period undulator has moderate gain without senous degradation due to 
beam quality. The CEBAF continuous beam of micropulses must only be 
slightly above the resonator loss to achieve saturation in strong optical fields. 
At normal incidence from 4 = 1000 -» 2000 A, clean Al in ultra-high vacuum has 
only 4% loss, and dielectric multilayers of MgF. on Al has only 2% loss [22]. 
The following three design modifications may improve gain, but with some 
added risk depending on the details of increased sensitivity to the electron 
beam quality. An increased sensitivity to beam quality may have greater than 
anticipated effect on gain because, until the accelerator is completed and 


begins operation, the true electron beam quality is unknown. 


Longer Undulator: When an FEL is not seriously degraded by beam 
quality, the design may benefit by increasing the undulator length L. As the UV 
undulator length is increased to L = 600 cm with N = 100 periods, the optimum 
Rayleigh length, Z,) = L/V12, can be used because the optical mode waist radius 
increases « L"*. The filling factor is then decreased to a more optimal value of 
F =0.5 with fewer electrons in the tails of the wider Gaussian optical mode. 
The dimensionless current density, j/F « N°F, is dramatically increased to 
jF = 11 and increases the potential for much higher gain. However, the penalty 
for a longer undulator is increased sensitivity to beam quality. The increased 
spreads due to energy spread and emittance are og=2.5 and o,=0.6. 
Assuming the energy spread is distributed as a Gaussian and the emittance is 
spread in the exponential distribution [16], the UV gain is calculated to be more 


than 70% for the longer undulator. This is an attractive gain for a FEL oscillator 
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at shorter wavelengths because mirror losses can be significant. With a longer 
undulator, it appears than the CEBAF UV FEL could be extended to reach 
wavelengths shorter than A =2000 A, but would become sensitive to the 


unknown shape of the f(v,) distribution. 


Smaller Electron Beam: The UV FEL gain is significantly reduced 
because some of the electrons are in the tails of the optical wavefront. Making 
the electron beam smaller in the transport system of the accelerator reduces 
this effect at the expense of an increase in the contribution of angular spread to 
poor beam quality in o,. If the radius of the beam at the waist is reduced to 
f, = 0.17 mm, then the angular spread is increased so that the beam radius 
roughly doubles its size to r, =0.3 mm at each end of the undulator. The 
previously negligible angular spread is now o,=1.2, and has a significant 
contribution on beam quality roughly equal to the energy spread contnbution, 
Og =1.3. The net result is favorable and raises the gain to 36% including 
diffraction, energy spread, and angular spread. As with the longer undulator, 
the smaller electron beam would give more gain, but would again become more 


sensitive to the unknown shape of the f (vo) distribution. 


FEL Klystron: The FEL klystron increases gain in weak optical fields at 
the expense of low-power saturation and increased sensitivity to beam quality 
[16]. The FEL klystron consists of two undulator sections separated by either a 
drift or a dispersive section. The drift and dispersive mechanisms are the same 
mathematically, but in practice, it is the dispersive section that gives a 
substantial increase to the FEL gain for a given interaction length. It is also 
advantageous that the strength of the dispersive section can be controlled 
dunng the experiment. Using the CEBAF UV FEL parameters including energy 


spread, the strength of the dispersive section was increased to find the 
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maximum gain. Because of the increased sensitivity to beam quality, no 
significant increase in gain was obtained over the conventional design. For a 
stronger klystron the effects of beam quality decreased the gain and countered 


the benefits of the klystron. 
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V. FELS WITH LARGE FILLING FACTOR 


A. INTRODUCTION 


As FEL experiments strive for shorter wavelengths, the size of the optical 
mode will become comparable to the size of the electron beam and the filling 
factor will increase. This was seen in the CEBAF UV FEL in Chapter IV. The 
theoretical single-mode gain is G = 0.135/F, but does not include the effects of 
self-consistent gain, beam quality, beam size or optical mode diffraction. 
According to the theoretical single-mode gain, as F increases, G increases. In 
reality, as the optical mode size approaches the electron beam size, fewer of 
the electrons are able to effectively participate in the bunching process. The 
result is that gain is reduced. Additionally, electrons across the electron 


distribution within the beam contribute to a larger phase shift. 


Figure 5-1 shows the size relation between a parabolic electron distribution 
and a Gaussian optical mode for a small, medium, and large optical mode. In 
the figure at the left, the optical mode is large and the electron beam is 
completely contained in the optical mode where the optical mode amplitude is 
large. All of the electrons will participate in the gain process. In the center 
figure, the optical mode and electron beam are of roughly equal size. Some of 
the electrons at the edge of the beam exist where the optical amplitude is small. 
These few electrons will not contribute significantly to the gain process. In the 
figure at the right, the electron beam is larger than the optical mode. Here, 
many of the electrons see a small optical amplitude and will not participate in 


the gain process. 
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Figure 5-1: Comparison of a parabolic electron distributicn (dark) to a large, 
medium and small Gaussian optical mode (light). 


B. THE PHASE VELOCITY RESONANCE SHIFT 
1. The Resonant Phase Velocity 


The maximum of the FEL gain spectrum is shifted by optical diffraction 
when Zp is finite [23]. As in the CEBAF UV FEL, a large electron beam size 


Causes an additional change in the phase velocity. 


To estimate the shift in resonance, consider the free Gaussian optical 
mode [16] as an approximate for the optical field when the gain ts low. The 


Gaussian optical field is 
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where r is the dimensionless distance off axis, w*(t) = 1+(t-ty)*/zé is the 
dimensionless optical mode radius, Z, is the dimensionless Rayleigh length, t is 
the dimensionless time along the undulator, t, is the position of the optical 
mode waist along the undulator, and a, is the dimensionless optical amplitude 
at r=0. From equation 5-1, a Rayleigh length of z, = 1/Vi2 will minimize the 


optical mode volume averaged over the length of the undulator. 


The electron pendulum equation [16, 21] describes the bunching of 
electrons in the presence of the optical field in equations 5-1. If the Rayleigh 
length is large, Z) > », equation 5-1 becomes a plane wave so that the optical 
mode is much wider than the electron beam, and the phase shift in equation 5- 
2 becomes negligible. When the Rayleigh length is finite and the electron beam 
is wide, the optical phase shift can be significant and can also influence the 
FEL interaction. To evaluate the effect of a finite Rayleigh length, equations 5-1 
and 5-2 are substituted into the pendulum equation ia lag (t)lcos[+d ¢(1)] 
where C = (k+ko)z(t}-wt is the electron phase derived in Chapter Ill. Expanding 
in powers of zo' to simplify the expression, to second order in zo', 
g(t) = tt, /Zo + *(t-1,,)*/Z6 + O(z5?) + +--+, the pendulum equation for an 
electron in a Gaussian optical mode with weak optical field strength, ao < zx, 
becomes 


: — | — 5-3 
Zo 226 Zo} Zo 
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where the electron phase is Cy = (k+ko)Z,—wt at position z and time t, k = 2x/2 is 
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the optical wavenumber, ky = 27/Ay is the undulator wavenumber, and G is the 


initial electron phase. The last term in equation 5-3 is 


eee 5-4 


where vo is the initial phase velocity. 


The amplitude of the field in equation 5-3 increases and decreases in 
time, t, along the undulator with maximum amplitude at the mode waist position 
t=t,. When the electrons are off-axis, the field amplitude is weaker. The 
electron phase is shifted by a term proportional to ty, but this shift is 
inconsequential because the electrons are randomly spread in €. The electron 
phase velocity is proportional to the evolving time +t, and is also shifted by the 


finite Rayleigh length and the radial position r. 


Equation 5-4 can be identified as the new phase velocity which has 
important consequences because maximum FEL gain occurs when Q = 2.6 [16]. 
For an FEL with a filament electron beam (r = 0 for all electrons), the maximum 
gain occurs at v™*=2.6 + 1/zZo. If 2) =1/¥12, the maximum gain occurs at 
Vo = 6 instead of vo>=2.6 for plane waves with Zz). Because the gain 


spectrum bandwidth is Avo = x, this is a significant shift [21]. 


lf the electron beam is not a filament, the effect of the electrons 
injected off-axis must be accounted for by finding the weighted average of the 
square of the off-axis injection distance over some normalized electron 


distribution function f(r). The distribution function is normalized such that 


ff (r) rdrd@ =1. The resonant phase velocity in equation 5-4 becomes 


ver n254 1) : 
0 


2 


Z0 


45 


where, in general, the average of the off-axis injection distance square is 
2x R 
<r> = J, J, r? f(r) rdrd®, where R is the maximum radius of the distribution. 
2. The Parabolic Electron Distribution 


If the electron beam cross-section has a parabolic shape so that 


fry=— 1 - | | 5-6 
76 @) 


2 2 
e 206 








with normalized width o, =/,(7/LA)* and maximum radius R = V20,, then the 


average square radius Is 


en Pee ee nie cee 
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Substituting equation 5-7 into equation 5-5, the resonant phase velocity for a 


parabolic electron beam of width o, is 


2 2 
vinx 964 1 _ a 


5-8 
Zo 326 





Simulations of the CEBAF UV FEL in Chapter IV showed that peak 
gain occurred at 2) = 0.45 and v» = 4.3. Equation 5-8 with z) = 0.5 and o, = 0.78 
predicts that maximum gain will occur at vo = 3.0, neglecting the effects of beam 
quality. For a Gaussian distribution, the phase velocity for peak gain is shifted 
by about og =4nNAy/y [16]. For the CEBAF case, og = 1.3. Including the effects 
of diffraction, beam size, and beam quality, the predicted value of 


Vo ~=2.6+2-1.6+ 1.3 = 4.3 in close agreement with the simulations. 


Figure 5-2 displays a graph of vg“ versus Rayleigh length for an FEL 
with low gain (j = 0.1) and a filament electron beam (o, =0.1). The dark line is 


a plot of equation 5-8, and the grey line is a plot of vg” including diffraction 
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only. The "X"s are data points taken from three-dimensional self-consistent 
simulations. The graph shows that for a small electron beam equation 5-8 is 
valid over a large range of Rayleigh lengths. The data points correspond 
closely with both plotted lines because the correction factor for the electron 
beam size, -202/z%, is small. The correction factor for the beam size is small 
as long as o, is small and zp is large. For small electron beams, equation 5-8 
is valid for z) > 0.2. Even where equation 5-8 is not an accurate estimate of the 
resonant phase velocity, the actual phase velocity tends toward the value 


predicted by equation 5-8. 


= 2.6+1/Z0-20¢7/3Z02% 
2.64+1/2Z20 


simulation output 
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Figure 5-2: Optimum phase velocity versus Rayleigh length for a small parabolic 
electron beam. Equation 5-8 is valid for z, > 0.2. 


Figure 5-3 displays a graph of vg’*™ versus Rayleigh length for an FEL 
with low gain (j = 0.1) and a large electron beam (o, = 0.7). Again the grey line 
iS a plot of vo™ including diffraction only. vg including the effects of diffraction 
and beam size is plotted as a dark line, and the "X"s indicate data points taken 


from simulations. This graph shows that for a large electron beam the 
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correction factor is larger and diverges from the actual vo™ at a larger Rayleigh 
length. The correction factor in equation 5-8 is larger so that equation 5-8 is 
valid over a smaller range of Rayleigh lengths. The difference can be attributed 


to the finite value of z, in the expansion for equation 5-3. 
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Figure 5-3: Optimum phase velocity versus Rayleigh length for a large 
parabolic electron beam. Equation 5-8 is valid for z, > 0.5. 


3. The Gaussian Electron Distribution 
Assuming a Gaussian distribution to the electron beam profile, then 


- 2 
= e r*r2a08 ; 5-9 


f(r) = ono? 


describes the electron distribution. The average electron squared radius is of 


_ 41 = -r/202 ne Se: 
eee J,re rdrd@=202 . 5-10 


Here the radius of the Gaussian distribution is taken to be infinite in order to 


complete the integration in equation 5-10. However, when performing 
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simulations using a Gaussian electron distribution a cutoff radius must be 
assumed. This will increase the errors of the final solution. Equation 5-10 is 
substituted into equation 5-3 to find the phase velocity for maximum gain 


including beam quality is 


2 

@) 
yeep =e), 5-11 
Zo Ze 


Overall, the equation for the resonant phase velocity including beam 
size is much less accurate for the Gaussian beam than for the parabolic beam. 
The reason is that many more of the electrons lie farther from the beam center. 
These electrons lead to a large value for the average square radius, but 
because the optical field is very weak at the edges, they contnbute little to the 
FEL gain and phase shift. The principle that a wide electron beam will lead to a 


shift in resonant phase velocity is still validated by this analysis. 


C. GAIN DEGRADATION WITH LARGE ELECTRON BEAMS 


Figure 5-4 shows the decrease in gain G(o,)/Gp as the electron beam 
increases beyond the optical mode waist [20]. The gain is determined 
numerically for a parabolic current density, and is compared to the theoretical 
single-mode gain Gy = 0.135/F for a filament beam, while o, is compared to wo. 
The Rayleigh length is z) = 0.5 and waist position is t, = 0.5 with phase velocity 
for near maximum gain at vo =4 in a weak optical field aj =1. The current 
density is taken to be j/F = 1.0. At each value of o,/wo the filling factor changes, 
and the peak current density / increases keeping jF fixed. At o,/w,. = 0.8, the 
actual gain, G(o,), and filament gain, Go, are equal. For larger beams, the gain 


decreases significantly below the filament gain. 
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Figure 5-4: Relative Gain versus Relative Beam/Mode size for low gain FEL. 


For values of jF < 1 down to jF =0.001, the curve in Figure 5-3 remains 


essentially unchanged. For jF 2 1, there is a small upward shift in gain on the 


left side of the curve at small values of o,/wo. For jF 2 10, the upward shift is 


more pronounced and is caused by optical mode distortion when the current j 


becomes large. For jF < 1, the value of z) can be varied from = 0.3 to =0.7 


with no significant change to the curve. The actual value of the current density 


for the CEBAF UV FEL is jF = 2.4 where j =2.7 and F = 0.9, so that the curve 


in Figure 5-4 can be used to describe the gain degradation in that experiment. 


Figure 5-4 appears to have a wide range of validity, and is useful for 


experimental design of FELs with large electron beams. 


50 


D. DISCUSSION 


The resonant phase velocity of an FEL is shifted by diffraction and by a 
finite electron beam size. By calculating the average off-axis electron injection 
distance, the new resonant phase velocity can be predicted over a range of 
Rayleigh ranges and electron beam sizes. The range of z. over which 
equations 5-8 and 5-11 hold is a result of the assumption that zp»1 in the 
expansion of the optical mode and the phase of the Gaussian optical mode. 
Because of this assumption, the analytical solution for vg“ is most accurate 
when 2, > 0.5. When longer Rayleigh ranges are used, the solutions above can 
be a sufficiently accurate estimate of the resonant phase velocity. For small 
electron beams (co, <0.5), the predicted resonant phase velocity remains 
accurate to smaller values of Z>. However, if the electron beam is large in 
comparison to the optical mode size, the solution for the resonant phase 
velocity at short Reyleigh ranges is less accurate, particularly for Gaussian 
electron bears. An empirical solution for the phase velocity for a parabolic 
beam is vo'™ = 2.6+1/z,-202/5z. The analogous empirical solution for 
Gaussian beam is v$™ = 2.6+1/z,-o2/zé. Both of these solutions are much 
more accurate over a larger range of electron beam sizes and Rayleigh lengths. 
They may be used to predict the resonant phase velocity for maximum gain in 


many experiments. 


The reduction in gain for an FEL with large filling factor results from the 
non-participation of electrons which lie near the edges of the electron beam. 
Figure 5-4 presented in Section C of this chapter, can serve as a kind of 
“universal” gain curve to provide a simple prediction for FEL gain for 
experiments with a wide range of dimensionless electron currents and filling 


factors. 
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Vi. OPTICAL MODE DISTORTION 


A. INTRODUCTION 


The optical field in an FEL can be thought of as the superposition of the 
input optical field and the excited optical field, which results from stimulated 
emission [24]. Assuming that the input optical field is a Gaussian, the field can 
be observed to focus toward the optical mode waist and then diffract out again 
away from the axis at a rate inversely proportional to the Rayleigh length. In 
dimensionless units Wo = YZo, and the optical mode radius is w* = 1+(1-1,,)/zé. 
‘The "new light” of the excited optical field begins near the electron beam. If the 
electron beam is much narrower than the optical mode, the excited field has an 
effective mode waist and Rayleigh length much smaller than the input optical 


mode and will diffract much more rapidly. 


B. THE NUMERICAL SIMULATION 


The effects of gain and electron beam size on the excited optical field can 
be observed through numerical simulations of the FEL interaction by subtracting 
the input optical field from the total optical field as they develop along the 
undulator. This chapter examines these effects with a modification of the 
simulation used to produce the output figures for the CEBAF UV FEL in 
Chapter IV. The simulation program uses numerical integration of the parabolic 
wave equation coupled with the electron Lorentz force equation [16]. This 
numerical simulation accounts for many transverse modes self-consistently and 


is general enough to include different undulator designs, optical mirror 


a2 


arrangements and driving currents. The evolution of the optical fields and 
electron currents are shown as an intensity/contour plot with the x and y axes 
scaled to VLA/x. The maximum field is white and the minimum field is black 


with two constant amplitude contours. 


These simulations begin in weak optical fields a) = 1 with phase curvature 
that gives Rayleigh length z)=0.3 that will focus at t, =0.5 without the FEL 
interaction. The electron beam will have a parabolic shape with radius, o,, 
phase velocity, v.=6, and current density, /, at the center of the beam. For 
these simulations it is assumed that there is no energy spread or emittance in 
order to observe the most simp > result of the optical mode distortion. Energy 


spread and emittance can be included when desired. 


The simulation follows two independent optical mode evolutions through 
the undulator. The first is the undistorted optical field developed with j =0. 
The second is the total optical field developed through the fully self-consistent 
simulation with j and o, equal to some appropriate value. The output of this 
simulation shows the difference between the total optical field amplitude and 
the input field amplitude as they develop along the undulator, Ala(x,t)!, and at 
the end of the undulator, Ala(x,y)l. Because this is a difference of amplitudes, 
it may display negative values. These negative values occur where the input 
optical field and the excited optical field destructively interfere so that the total 
optical field is less than the input optical field. The difference in phase between 


the total optical field and the input optical field along the axis is plotted as A(t). 
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C. THE INPUT OPTICAL FIELD 


Figure 6-1 shows a Gaussian optical mode in an FEL undulator. The 
dimensionless current density, /, is zero so there is no gain, and the sever 
(lower-right) remains constant. Because there is no excited field, the input 
optical mode remains undistorted throughout the undulator. In the evolution of 
the optical field amplitude through the undulator, |a(x,t)! (upper-left), the optical 
mode focuses at the optical mode waist, ty, and later begins to spread. 
Because the dimensionless Rayleigh length is z)=0.3, the optical mode area 
doubles at At=+0.3 on either side of t,. The input optical mode amplitude at 
the end of the undulator, la(x,y)l (upper-middle), remains in a Gaussian 
distribution. The blocks which describe the electrons, o(x,t) and o(x,y), and the 
phase velocity evolution, f(v,t), are blank because there are no electrons used 


in the simulation. 


D. THE EFFECT OF BEAM SIZE ON MODE DISTORTION 


The waist of the excited optical field, w,, is roughly equal to the electron 
beam size [24], so that the Rayleigh length of the new light can be much 
shorter than that of the input optical mode. Therefore, the Rayleigh length for 
the new light is determined by z, = o%. If the input optical mode is much larger 
than the electron beam then z,) > z,, and the excited optical field will diffract 
much more quickly than the input optical field. To show the effect of the 
electron beam size on the excited mode, several simulations were run where 
the product jo? is held constant in order to maintain the same single mode 


gain. Figures 6-2 through 6-5 show four examples where joé = 0.1. 
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Figure 6-1: Free Gaussian optical mode in an FEL undulator. 
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In Figure 6-2, the wide electron beam results in a wide excited mode, 
Ala(x,t)!|. The excited mode remains close to the electron beam axis. The 
difference in the phase between the total optical field and the input optical field, 
Ad(0,t), is very small. In Figure 6-3 the excited optical field is initially more 
narrow, but diffracts to a larger area at the end of the undulator. The area just 
Outside the second contour corresponds to nearly zero amplitude. Farther 
outside the second contour there is destructive interference between the input 
and excited modes so that the resulting difference in amplitudes is negative. 
The phase difference between the two modes has grown, but is still small. 
Figure 6-4 continues the trend as the waist continues to narrow causing rapid 


diffraction with a small phase difference. 


In Figure 6-5 there is a major difference. The overall trend of the 
narrower electron beam rapidly diffracting away from the undulator axis is 
continuing, but a hole has developed at the center of the excited optical field 
* near the end of the undulator. The size of the hole roughly corresponds to the 
size of the electron beam. The evolution of the average electron dnving 
Current, o(x,t) shows that the electrons are bunched at the end of the undulator. 
The gain evolution, In(i+G(t)), shows a maximum at the end of the undulator. 
The excited field at the end of the undulator in Figure 6-5 forms a "doughnut" 


around the undulator axis. 
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i sg 6-2: The excited optical field of a low gain FEL with a wide electron 
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Figure 6-3: The excited optical field of a low gain FEL with optical mode slightly 
larger than the electron beam. 
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Figure 6-5: The excited optical field of a low gain FEL with a small electron 
beam. A hole has developed in the excited optical field. 
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The explanation of this hole involves the relative phase between the new 
light emitted along the undulator. In the development of equation 5-3 in 


Chapter V, the phase of the Gaussian optical mode was expanded to 


7 2 
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At the undulator axis (r = 0) the relative phase between the optical mode at any 


two points separated by At is 


Adc = st . 6-2 


0 


The effective mode waist for the excited field is roughly equal to the 
electron beam size so that, z)=0o7%. The phase shift in equation is then 
Adg =Avoé. When Adg =x the excited optical field is eliminated by destructive 
interference when o, = V(Awn). The first light emitted near t = 0.6 destructively 


interferes with light emitted at the end of the undulator, t= 1, when o, = 9.36. 


Figures 6-6 and 6-7 show two additional examples of this destructive 
interference phenomenon. These figures, unlike Figures 6-1 through 6-7, are 
plotted with a greyscale which runs from white to black with increasing 
magnitude and no contours. This allows a better display of the thin, excited 
optical mode and electron beam. Also, the transverse window width of the 
displays is 8 instead of 4 because of the increased diffraction by the narrow 


electron beam. 


In Figure 6-6 the simulation is run with /=10 and o, =0.1, so that the 
product of j o2 remains 0.1 for low gain. The destructive interference is seen 


along the undulator axis extending along the final third of the undulator. 
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Figure 6-6: The excited optical field of a low gain FEL with electron beam size 


6, =.1. The area of destructive interference extends over the final third of the 
undulator. 
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Figure 6-7 is the same type of output with / = 40 and o, =0.05 so that 
jo2=0.1 again. Here, the area of total destructive interference extends through 
the entire second half of the undulator. This shows that the simple analysis 


above is correct. 


The above explanation of the destructive interference is simple but valid. 
Actually, light is continuously emitted along the undulator and the amplitude of 
the light emitted is proportional to the gain. A complex integral is required to 
fully explain the phase relation that leads to destructive interference. However, 
the simple argument shows that as the electron beam radius decreases the 
area of destructive interference within the excited optical mode will extend 


toward the front of the undulator. 


E. THE EFFECT OF GAIN ON MODE DISTORTION 


In order to examine the effect of gain on optical mode distortion, the 
electron beam size is held constant while the dimensionless current density is 


increased. Figures 6-8 and 6-9 are two examples. 


Figure 6-8 shows a low gain system with / =0.1 and o, =0.5. The final 
gain is only = 0.5% and the power remains roughly constant over the length of 
the undulator. The maximum amplitude difference between the input optical 
mode and the total optical mode is small (Ala|,,,, = 0.002), as expected from low 
gain. With the moderate electron beam size, the excited optical amplitude 
difference diffracts with a Rayleigh length z, = 1/V2. This is roughly twice the 
input optical Rayleigh length z)>=0.3 so the excited mode spreads at roughly 


half the rate of the input mode. The phase difference along the axis is small. 


63 


See 


rn NS 
EN 


oo 


4 ? oP Sete ‘5 A a tee “ELE ie 
v ay SR RAS q ¢ 5 ETE NS ot GULLS OPIS IIILIESTTD 
Ce ee ae 


SY 


Cebit senna 


> 5 


-. a 
PERI. 
Ves ERS ‘ (ict epees 
BE 
PRE Zi SEO TE GIES: Gp, 
a igh Fs i cf 
() 


le 


tt) 


LE MAZE 


Aa(x,t) 0.025 
-0.018  Aa(x,y) 0.026 
0 G(x,t) 1.8x104 
0 o(x,y) 1.8x104 








Ad (0,T) 


Ln (1+G(T) ) 


Figure 6-7: The excited optical field of a low gain FEL with electron beam size 
o, = .05. The area of destructive interference extends over the entire second 
half of the undulator. 
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Figure 6-9 shows the excited mode of a high gain system. Here the 
electron beam is the same as in Figure 6-8, but the dimensionless current 
density is three orders of magnitude larger. At the lower-nght of Figure 6-9, the 
final gain is 18.5% and the power has grown significantly. The profile of the 
excited mode at the upper left again shows diffraction at roughly half the rate of 
the input mode. 

It is interesting to note that the plots of the excited optical modes in 
Figures 6-8 and 6-9 apper to be almost identical. The contours of Figure 6-8 
appear somewhat longer, but the shapes are similar. There is a difference in 
the optical phase evolution in thy two plots. In Figure 6-8, the phase difference 
is very small (Ad = 0.0046), so there is nearly maximum constructive interference 
along the undulator axis. In Figure 6-9 the phase difference is still small 
(Ad = 0.5), but more significant. The interference between the input optical 


mode and the excited optical mode along the undulator axis is not maximized. 


F. DISCUSSION 


The three-dimensional simulation displaying the excited optical field 
amplitude and phase provides a useful tool for examining optical mode 
distortion. The effective Rayleigh length of the excited optical field is = o% and 
the diffraction of the excited optical field is dependent upon the beam radius. 
The destructive interference caused by the phase difference, Ao = Avo; = x, 
when At= 1/3 occurs when o, = 0.3 and leads to the development of a hole 


within the excited optical mode. 
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| 0 O(x,T) 3.4 
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Figure 6-8: The excited optical field of a low gain FEL. 
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Figure 6-9: The excited optical field of a high gain FEL. 
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With constant electron beam size, the shape of the excited optical field 
appears to be relatively independent of the FEL gain. For both a low gain 
system (j=0.1) and a high gain system (j = 100) the shape of the excited 


optical mode is roughly constant. 


Further use of this tool to study the excited optical mode may lead to 


additional, more significant, insights into optical mode distortion. 
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Vil. CONCLUSIONS 


The development of FELs as a source of coherent radiation is continuing 
as several theoretical and technological issues are met and resolved. The 
military applications of the FEL, whether on the ground, in space, or aboard a 
ship has several significant advantages over conventional kinetic systems and 
other high-energy laser systems. The rapid response, wavelength tuneability, 
and "infinite" magazine make the FEL a highly desirable system for shipboard 
use. The all-electric operation of the FEL is especially attractive when coupled 
with the design of the Navy’s all-electric ship. With these important attributes, 
research should continue toward the development of high power FELs with the 


goal of shipboard application. 


The proposed CEBAF UV FEL represents a significant demonstration of 
the technologies required for a high power FEL. The superconducting 
accelerator, high-average power, and short wavelength of the CEBAF UV FEL 
will reach the highest levels of performance ever demonstrated. Simulations of 
the proposed experiment presented in Chapter IV helped in the FEL design 
selection. The designs presented appear feasible and should result in 
successful accomplishment of the design goals. 

In short-wavelength FELs, such as the CEBAF UV FEL, the optical mode 
is small and can be roughly the same size as the electron beam. The shift in 
the resonant phase velocity resulting from the injection of electrons away from 
the undulator axis is predicted in Chapter V. In some instances this shift can 


be significant. There is also an inevitable reduction in gain resulting from the 
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off-axis electrons interacting with the weaker optical fields. The gain of sucha 
system can be predicted with the aid of the universal curve also presented in 


Chapter V. 


Optical mode distortion is a key topic in the understanding of the FEL. A 
new tool presented in Chapter VI helps examine and understand this 
phenomenon. The modified three-dimensional FEL simulation displays the 
excited optical mode in the undulator. The size of the electron beam has a 
significant effect on the charactenstics of the excited optical mode. For very 
small electron beams the phase difference between new light emitted near the 
axis along the undulator can result in total destructive interference within the 
excited optical mode. This interference can cause a “dimple” to occur at the 
center of the total optical field. This dimple and other effects in the excited 
optical mode could be researched further with the aid of the modified 


simulation. 
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